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The quark-gluon vertex in Landau gauge QCD is investigated in the 
Dyson-Schwinger approach. The aim is to obtain a fully self-consistent so- 
lution of the quark propagator and the quark-gluon vertex Dyson-Schwinger 
equation (DSE) in the vacuum using the gluon propagator as input from 
other calculations. The truncation scheme used in an earlier study is sys- 
tematically improved and the first decisive step towards a full solution is 
presented. First results for the quark propagator with a fully dressed quark- 
gluon vertex in a truncation that incorporates all twelve vertex tensor- 
structures are shown. 

PACS numbers: 11.15.-q, 11.30.Rd, 12.38.Aw 

1. Motivation 

Prominent features of QCD like confinement and dynamical chiral sym- 
metry breaking (D^SB) are still not satisfactorily understood. One possi- 
bility to explore the non-perturbative regime of QCD is to employ DSEs, an 
infinite tower of coupled integral equations, see e.g., ref. pQ and references 
therein. Their non-perturbative nature makes them an appropriate tool for 
exploring the low-energy domain of the theory. Careful truncations have 
to be imposed on the system in order to evaluate the equations. In these 
proceedings we will extend a prior study of the quark-gluon vertex [2] by 
systematically improving on the used truncation. The aim is to shed more 
light on a possible relation between confinement and D^SB by studying the 
quark-gluon vertex. 

The starting point is the system of equations for the quark propagator 
and the quark-gluon vertex in Landau gauge using fits for the gluon propa- 
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gator as input [21 [3] . The three-gluon vertex is modeled as in ref. [2jf]The 
infrared suppression of the gluon propagator relates on the one hand to an 
infrared enhancement of the ghost propagator (see e.g., ref. [6]) and on the 
other hand to the positivity violation for transverse gluons. The latter prop- 
erty indicates the role of the transverse gluons in a non-perturbative BRST 
quartet mechanism [7] and therefore a certain aspect of gluon confinement: 
In Landau gauge QCD transverse gluons are confined and therefore not 
confining. Thus, given the infrared suppression of the gluon propagator, 
the quark-gluon vertex has to acquire a certain strength in the infrared to 
provide D^SB and, eventually, quark confinement. In ref. [2] it has been 
found that within the imposed truncation (and at vanishing temperatures 
and densities) D^SB and quark confinement occur either together in a self- 
consistent way or are both absent. By improving on the truncation of the 
quark-gluon vertex DSE we hope to get deeper insights into the possible re- 
lation between quark confinement and D^SB. In addition, isolating impor- 
tant tensor structures of the quark-gluon vertex will likely lead to improved 
phenomenological applications. Especially, corresponding studies at non- 
vanishing temperatures and/or quark chemical potentials can benefit from 
a more detailed understanding of the quark-gluon vertex, see e.g., ref. [8]. 



2. The coupled system 

The coupled system of the quark propagator and the quark-gluon vertex 
is shown in Fig. [T] This system serves as a starting point for this study. The 
vertex equation in Fig. [l]has been derived from a 3PI effective action and 
corresponds to the DSE truncated at three-point level [9J. Note that as usual 
in an nPI-based approach all vertices are dressed. For obvious reasons, the 
last two diagrams on the right-hand side of the vertex equation are referred 
to as the non-Abelian and the Abelian diagram. Upon performing the color 
traces one finds that the non-Abelian diagram has a prefactor of N c /2, 
whereas the Abelian diagram is suppressed by a factor of — l/2iV c . This 
suppression has been found to be even more severe due to the dynamics 
of the system [;2j. Thus, as a first step to solve the system of the quark 
propagator and the quark-gluon vertex DSE, only the non-Abelian diagram 
has been taken into account. 



1 These fits relate to the so-called scaling solution, one of two types of solutions for 
Landau gauge Green functions when classified by their infrared behavior, see e.g., 
refs. [HE] an< i references therein. However, we want to emphasize here the evidence 
that the difference between these types of solutions is irrelevant for phenomenology. 
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Fig. 1. The quark propagator and the quark-gluon vertex DSE. All internal prop- 
agators are dressed. The vertex equation has been derived from a 3PI effective 
action [§]. 

2.1. A first step 

The system as depicted in Figs.[2]and[3]has been considered as a first step 
towards the full solution of the system of Fig. [TJ Within this approximation, 
one vertex of the non-Abelian diagram remains bare while the second is 




Fig. 2. Quark propagator DSE. Full blobs denote fully dressed quantities, shaded 
blobs arc fits from other studies. 




Fig. 3. The truncated quark-gluon vertex DSE. Full blobs denote fully dressed 
quantities, shaded blobs are fits from other studies. For the three-gluon vertex a 
model has been used which is adopted from [2]. 



dressed such that the full vertex experiences feed-back from all twelve tensor 
structures. This vertex also enters the propagator equation Fig. [2j In the 
vacuum, twelve tensor structures are needed as a basis to span the vertex. 
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Instead of employing the widely-used Ball-Chiu basis [10] it turned out to 
be numerically advantageous to use a simpler basis given by 



r M (pi,P2) oc 



1 

Pi 




where p\ and p 2 are the in- and outgoing quark momenta. As input from 
the Yang-Mills sector the gluon propagator has been taken from earlier 
DSE studies^] The three-gluon vertex is represented by its tree-level tensor 
structure combined with the following ansatz [2] 
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where x = pi + p\ + Vz is the sum over the squared gluon momenta entering 
the vertex and k ~ 0.595. This ansatz reflects the strong IR enhancement 
known from Yang-Mills theory (see, however, the remark in sect. 3) as well 
as ensures the correct running in the UV region. The residual parameters d\, 
d 2 and ds have to be fixed to physical observables, e.g. the chiral condensate. 



2.2. Obtaining the vertex dressing functions 

As a first step, the dressing functions have to be disentangled to get them 
in an explicit form. Since the employed basis is neither orthonormal nor 
orthogonal one obtains a linear system of equations for the vertex dressing 
functions in terms of the twelve right-hand side projections. This system 
can be solved in advance, giving rise to explicit expressions for the dressing 
functions as linear combinations of the right-hand side projections^} In 
each iteration step these projections are evaluated numerically and are used 
subsequently to obtain the dressing functions by inserting them into the 
solution of the pre-calculated linear system of equations. 



2.3. Renormalization and numerical treatment 

The quark wave function renormalization constant the mass renor- 
malization constant Z m as well as the quark-gluon vertex renormalization 
constant Z\p has been fixed within a MOM scheme. In particular, the 
renormalization conditions for the quark dressing function A{fi 2 ) = 1 and 

2 Here, a scaling-type fit obtained from a self-consistent solution of the corresponding 
ghost-gluon system has been used [3J. The fit parameter values are taken from [2]. 

3 The algebraic manipulations have been carried out with the program FORM 
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the tree-level vertex dressing function Ai(// 2 , /i 2 , 2/i 2 ) = 1 have been em- 
ployed, where fi 2 = 170 GeV 2 denotes the renormalization scale. Only the 
tree-level tensor 7^ has to be concerned in this procedure since all other 
tensor structures lead to UV finite expressions. Z m has been fixed by the 
condition B([i 2 ) = thr, where rriR is the mass at the renormalization scale. 

The evaluation of the coupled system of the quark propagator and the 
quark-gluon vertex as depicted in Fig. [2] and Fig. [3] is numerically de- 
manding. The actual calculations are thus performed on a cluster featuring 
Graphics Processing Units (GPUs). GPUs have already been applied suc- 
cessfully to solve DSEs |12j . They will be vital once the full non-Abelian 
diagram and/or the Abelian diagram is taken into account. 



3. First results 



The system described in Sec. 2.1 has been solved self-consistently, tak- 
ing a ll twelve tensor structures for the quark-gluon vertex into account. Fig. 
4(a) shows the behaviour of the mass function M(p 2 ) as obtained from the 
quark propagator dressing functions in the chiral limit. Due to stability 
issues it was not yet possible to ta ke the full IR strength of the three-gluon 
vertex (p 2 )~ 3ft into account]^} Fig. 4(b) shows the vertex dressing function 




Fig. 4. Left: The mass function M (p 2 ) — B(p 2 )/A(p 2 ) as obtained from the coupled 
system described in Sec. 2.1 Right: The vertex dressing function Ai (p 2 ) of the tree- 
level part 7^ of the quark-gluon vertex evaluated at p\ — p\ — p 2 and p\.p2 = 0. 



Ai which corresponds to the tree-level structure 7^ evaluated at the symmet- 
ric point p\ = p\ = p 2 and p\.p2 = 0. The significant IR enhancement seen 
in this function is also present in most other tensor structures, independent 
whether they have been generated by D^SB or are chirally symmetric. 



4 In our calculations the corresponding IR exponent has been varied from to —2k 
resulting in moderate changes of the dynamical generated mass. 
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4. Conclusions and outlook 

The coupled DSEs for the quark propagator and the quark-gluon vertex 
in Landau gauge has been investigated. Results for the quark propagator 
coupled to a partially dressed quark-gluon vertex including all twelve tensor 
structures have been presented. Isolating the important tensor structures 
will become vital for subsequent calculations taking the full non-Abelian 
as well as the Abelian diagram into account, aiming towards a complete 
solution of the system as depicted in Fig. [TJ Such a reduction of complexity 
is also mandatory in studies involving the quark-gluon vertex, especially 
those at non- vanishing temperatures and/or quark chemical potentials. 
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